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Abstract 

o ; 

OA , In this article, we study the strong decays of the newly observed charmed mesons 

J> ■ £>(2550), L>(2600), L>(2750) and L>(2760) with the heavy quark effective theory in the 

leading order approximation, and tentatively identify the (D(2550), D(2600)) as the 25 
doublet (0", 1") and the (D(2750), D(2760)) as the ID doublet (2", 3"), respectively 
The identification of the D(2750) and Z?(2760) as the same particle with J p = 3~ is 
disfavored. 



(N 



(N 
> 



PACS numbers: 13.25.Ft; 14.40.Lb 

Key Words: Charmed mesons, Strong decays 



1 Introduction 



Recently the Babar collaboration observed four excited charmed mesons -D(2550), D(2600), 
£>(2750) and£>(2760) in the decay channels D° (2550) -»■ D* + tt-, L>°(2600) -> D* + Tr~, D+vT 
£>°(2750) -> D* + ir-, L>°(2760) -»■ D+vr", L>+(2600) -> L>°7r+ and £>+(2760) D°vr + re- 
spectively in the inclusive e + e~ — > cc interactions at the SLAC PEP-II asymmetric-energy 
collider pQ, see Table 1. The Babar collaboration also analyzed the helicity distributions to 
determine the spin-parity, and suggested that the (D(2550), D(2600)) (denoted as (D', D*') 
respectively in Table 2) may be the 25 radial excitation of the (D,D*), and the .D(2750) 
and D(2760) may be the D-wave states. Furthermore, the Babar collaboration measured 
the following ratios of the branching fractions: 

X : Br (g(246o)° -> py) _ 147± 03 ±016j 
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: 0.32 ±0.02 ±0.09, 

= 0.42 ±0.05 ±0.11. (1) 

In the heavy quark limit itiq — > oo, the heavy- light mesons Qq can be classified in 
doublets according to the total angular momentum of the light degrees of freedom sg, 
H = Sq + L, where the Sq is the spin of the light antiquark q and the L is the orbital 
angular momentum of the light degrees of freedom [21 [3]. In the quark models, we usually 
use the n to denote the radial quantum number. In the case n = 1, for L = 0, the 
doublet (P,P*) have the spin- parity jf = (0~,l - )i; L = 1, the two doublets (Pq,P{) 

*- 2 

and (PijPo) have the spin-parity = (0 + ,l + )i and (l + ,2 + )3 respectively; L = 2, the 

2 2 

two doublets (P*,P2) and (P^j-fb) have the spin-parity = (l~,2~)s and (2~,3~)5 

*■ 2 2 
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Mass [MeV] 


Width [MeV] 


Decay channel 


L> u (2550) 


2539.4 ±4.5 ±6.8 


130 ± 12 ± 13 


D*+tt- 


L> u (2600) 


2608.7 ±2.4 ±2.5 


93 ± 6 ± 13 


D + tt-,D* + tt- 


L> u (2750) 


2752.4 ±1.7 ±2.7 


71±6±11 


D* + TT~ 


L> u (2760) 


2763.3 ±2.3 ±2.3 


60.9 ±5.1 ±3.6 


D + TT- 


L>+(2600) 


2621.3 ±3.7 ±4.2 


93 


D U 7T+ 


L>+(2760) 


2769.7 ±3.8 ±1.5 


60.9 


D U 7T+ 



Table 1: The experimental results from the Babar collaboration. 



respectively; where the superscript P denotes the parity. The n = 2, 3, 4, • • • states are 
clarified by the analogous doublets, for example, n = 2, L = 0, the doublet (P' , P*') have 
the spin-parity J? = (0 _ ,l-)i; n = 2, L = 1, the two doublets (P$',P{') and (Pi,P^) 

have the spin-parity jf t = (0 + ,l + )i and (l + ,2 + )3 respectively. 

The helicity distributions favor identifying the D°(2550) as the 0~ state, the Z)°(2600) 
as the 1~, 2 + , 3 _ state, and the L>°(2750) as the 1 + , 2~ state pp. From the Review of 
Particle Physics [4j, we can see that only six low-lying states, D, D*, Do (2400), Z>i(2430), 
Z?i(2420) and Z?2(2460) are established, while the 2S and ID states are still absent. 
The newly observed charmed mesons -D(2550), -D(2600), -D(2750) and -D(2760) may be 
tentatively identified as the missing 2S and ID states. 

The mass is a fundamental parameter in describing a hadron, in Table 2, we present 
the predictions from some theoretical models, such as the relativized quark model based on 
a universal one-gluon-exchange-plus-linear-confinement potential [5], the semirelativistic 
quark potential model [6], the relativistic quark model includes the leading order 1/Mh 
corrections [7], the QCD-motivated relativistic quark model based on the quasipotential 
approach [8], for comparison. From the Table, we can see that the masses of the -D(2550), 
£>(2600) and D(2750), L>(2760) lie in the regions of 2S and ID states, respectively. 

In Ref.[9], Sun et al study the strong decays of the D(2550), L>(2600) and L>(2760) in 
the 3 P model, and identify the L>(2600) mixture of the 2 3 5i — l 3 D\ states and the 
D(2760) as either the orthogonal partner of the -D(2600) or the l 3 .D3 state. In Ref.[10j. 
Zhong studies the strong decays of the D(2550), D(2600) and -D(2760) in a chiral quark 
model, and identifies the £)(2760) as the 1 3 -D3 state and the D(2600) as the low-mass 
mixing state of the 1^D\ — 2 3 Si states. 

In this work, we study the strong decays of the newly observed charmed mesons with 
the heavy quark effective theory in the leading order approximation to distinguish the 
different identifications. There have been several works using the heavy quark effective 
theory to identify the excited D s mesons, such as the D s (3040), D s (2700), D s (2860) 

puna nana. 

The article is arranged as follows: we study the strong decays of the newly observed 
charmed mesons with the heavy quark effective theory in Sect. 2; in Sect. 3, we present the 
numerical results and discussions; and Sect. 4 is reserved for our conclusions. 
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Experiment [TJ H] 


GI 15 


MMS [6 


PE [71 


EFG [8] 


D 


15^0" 


1867 


1880 


1869 


1868 


1871 


D* 




2008 


2040 


2011 


2005 


2010 


Dn 


ip Ao + 

2 


2400 


2400 


2283 


2377 


2406 


D[ 


lPil+ 


2427 


2490 


2421 


2490 


2469 


Di 


1P|1+ 


2420 


2440 


2425 


2417 


2426 


n* 

u 2 


1P§2+ 


2460 


2500 


2468 


2460 


2460 


D\ 




7 2763 


2820 


2762 


2795 


2788 


D 2 


1D|2- 


7 2752 




2800 


2833 


2850 


D> 2 * 


1D§2- 


7 2752 






2775 


2806 


D 3 


1D|3- 


7 2763 


2830 




2799 


2863 


D' 


2S i0" 


7 2539 


2580 




2589 


2581 


D*' 


25^1" 


7 2609 


2640 




2692 


2632 



Table 2: The masses of the charmed mesons from different quark models compared with 
experimental data, and the possible identifications of the newly observed charmed mesons. 



2 The strong decays with the heavy quark effective theory 

In the heavy quark effective theory, the spin doublets can be described by the effective 
super-fields H a , S a , T a , X a and Y a , respectively [15J, 
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where the heavy field operators contain a factor y/Mp and have dimension of mass |. The 
ground state and radial excited state heavy mesons with the same heavy flavor have the 
same spin, parity, time-reversal and charge conjunction properties except for the masses, 
and can be denoted by the super-fields: Ha, H a , H a , ■■■ ', S a , S a , S a , ■■■ ; 

T a , T' T" 

• • • ; etc, where the superscripts /, II and III denote the first, the second and the third 
radial excited states, respectively. With a simple replacement of the components P a , P*, 



P, 



Oct ■ 



to the corresponding radial excited states P a ', P* , PqJ, 



we can obtain the 



corresponding super-fields H' a , S' a , 
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The light pseudoscalar mesons are described by the fields £ = e fa , where 



A4 



2 ' 

7T~ 



At the leading order, the heavy meson chiral Lagrangians Ch, Cs, Ct, £-x, £-y for the 
strong decays to D^tt, D^-q and D { * ] K are written as [T6l fT7| [T8| [T9l [20] : 

£// = fi^Tr {#^7^75^} , 

£s = gs^ {HaSbj^A^} + h.c. , 

C T = ^Tr {BaT^iD^A + ipA^bajsj + h.c, 
A x 

C x = ^Tr{H a X^(iD^A + iPA)bal5} + h.c, 
A x 

C Y = -j^Tr {B a Y b ^[h{D^D v }Ax + ^{D^D x A v + D u D x A^^}+h.c.,(3) 

X 

where 

A* = , (4) 



A x is the chiral symmetry-breaking scale and taken as A x = 1 GeV [llj . the strong 
coupling constants gn, 9s, 9t, 9x and gy = (k\ + can be fitted phenomenologically if 
there are enough experimental data. The subscript indexes H, S, T, X and Y denote the 
interactions between the super-field H and the super-fields H , S, T, X and Y, respectively. 
We have smeared the superscripts /, //, ///, • • • for simplicity, the notation gu denotes 
the strong coupling constants in the vertexes HHA, H'HA, H'H'A, H"HA, • • • , the 
notations gs, gT, 9x and gy should be understood in the same way. In this article, we 
intend to study the ratios among different decay channels, the strong coupling constants 
are canceled out with each other, and cannot lead to confusion. 

From the heavy meson chiral Lagrangians Ch, £-s, £-t, £x, £-y, we can obtain the 
widths r for the strong decays to jDW-tt, D^rj and D^K easily, 

T _ Pern m 2 (r \ 

where the T denotes the scattering amplitudes, the p cm is the momentum of the final 
states in the center of mass coordinate. 

In calculations, we take the approximation A^ ~ i 9 y M . In the case that the light 
pseudoscalar meson momenta are not very small, we should add other terms and introduce 
new unknown coupling constants. Furthermore, the flavor and spin violation corrections 
of order 0(1/?h,q) to the heavy quark limit may be sizable, again we should introduce 
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new unknown coupling constants, which will not necessarily canceled out in the ratios of 
the decay widths. We cannot estimate the role and the size of such corrections on general 
grounds, however, we expect that they would not be larger than (or as large as) the leading 
order contributions. 



3 Numerical Results 

The input parameters are taken from the particle data group M 7r + = 139.57 MeV, M w o = 
134.9766 MeV, M K + = 493.677 MeV, M v = 547.853 MeV, M D + = 1869.60 MeV, M D o = 
1864.83MeV, M D + = 1968.47MeV, M D ,+ = 2010.25 MeV, M D *o = 2006.96 MeV, M D *+ = 
2112.3 MeV, M D( t, m = 2460.1 MeV @]. 

The numerical values for the widths of the strong decays 

D* 2 -> D* + ir~, D + tt- , 

D' -)• D* + ir~, D*°TT°, 

D*'(Dl, D 2 , D 3 ) -»• D* + tt~, D+TT-, D*°tt°, D°ir°, D*°r], D°r], D* s + K~, D+K~ ,(6) 

are presented in Tables 3-4. 

In Table 5, we present the experimental data for the ratio y-(d*+tt~) oi ^ ne we ^ es ~ 
tablished meson 1^(2460) from the Babar CLEO [ZE1[22], ARGUS [23], and ZEUS 
[24J collaborations, the prediction 2.30 from the heavy quark effective theory in the lead- 
ing order approximation is in excellent agreement with the average experimental value 
2.35. Compared with the experimental data from the Babar collaboration y(d*+tt~) = 
1.47 ± 0.03 ± 0.16 m, the heavy quark effective theory in the leading order approximation 
leads to a larger ratio. 

The total decay widths of the {D (2550), D (2600)) with the spin-parity (0~,l~)i are 

Td' ~ GeV and Yd*i ~ 2.0g^ GeV, the ratio ~ 0.85, which is smaller than 

the experimental data p- 2 ^- = 1.40, where we have used the central values of the widths 
r D / « (130 ± 12 ± 13) MeV and T D */ = (93 ± 6 ± 13) MeV from the Babar collaboration 
[T]. For the charmed mesons, the leading power flavor and spin violation corrections (of 
order 0{l/rriQ)) to the heavy quark limit may be sizable, we have to introduce new un- 
known coupling constants, the discrepancy may be smeared with the optimal parameters, 
furthermore, more precise measurements are needed to make a reliable comparison. In 

p 

the case of the ratio f—^, the prediction 0.30 from the heavy quark effective theory in 

D* 

the leading order approximation is also smaller than the experimental data 0.48 from the 
Review of Particle Physics [1], if the leading power spin corrections to the heavy quark 
limit are taken into account, the discrepancy can be smeared [25]. 

The ratio r(D*'->D*+Tr-) = *-*-82 from the heavy quark effective theory in the leading 
order approximation is larger than the experimental data 0.32 ±0.02 ±0.09 from the Babar 
collaboration [lj, just like in the case of the ratio r(D*^>D*+n~) ' anc ^ a § am more precise 
measurements are needed to make a reliable comparison. The strong coupling constants 
9d*Dtt arid Qd*d*-k receive sizable contributions from the flavor and spin violation correc- 
tions [201 [26], in the present case, the strong coupling constants gD*'D-n and gr>*'D*-K also 
receive the flavor and spin violation corrections besides the leading order strong coupling 
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constant gn, which maybe account for the discrepancy. We can tentatively identify the 
(£>(2550),£>(2600)) as the doublet (0~, l~)i with n = 2. 

2 

The existing theoretical estimations for the strong coupling constant gn among the 
ground state heavy mesons (n = 1) vary in a large range gn = 0.1 — 0.6, it is difficult to 
select the ideal value (one can consult Ref.[27] for more literatures), we usually use the 
value determined from the precise experimental data on the decay D* + — > D°7t + from the 
CLEO collaboration [28, 29J. In the present case, the strong coupling constants involve the 
radial excited S'-wave heavy mesons and ground state D-wave heavy mesons, therefor the 
situation is more involved, and it is impossible to determine the revelent parameters with 
the heavy quark effective theory itself without enough experimental data. The theoretical 
works focus on the strong coupling constants gu, <7S, 9t of the ground state S-wave and 
P-wave heavy mesons (one can consult Refs.pO j 127 } I30j for more literatures), while the 
works on the strong coupling constants gH, 9s i 9t of the radial excited S'-wave and P- 
wave heavy mesons and gx, gy of the ground state D-wave heavy mesons are rare due to 
lack experimental data |31| . In this article, we take the strong coupling constants gn, gr, 
gx and gy as unknown parameters, and prefer the ratios of the decay widths in different 
channels to compare with the experimental data. 

From Table 4, we can see that if we identify the (.D(2760), D(2750)) as the doublet 

(l~,2 - )3 with n = 1, the ratio r(D2->D*+n-) = fr° m the leading order heavy quark 
effective theory deviates from the experimental data 0.42 ± 0.05 ±0.11 greatly [iffi which 
requires the flavor and spin violation corrections depressed by the inverse heavy quark 
mass 1/rriQ are as large as the leading order contributions and have opposite sign, it is 
impossible, as the heavy quark effective theory has given many successful descriptions of 
the hadron properties [2 13 ED]. On the other hand, if we identify the (£>(2750), L>(2760)) 
as the doublet (2~,3~)s with n = 1, the deviation of the ratio t^w-t^t -\ = 0.80 
from the upper bound of the experimental data 0.42 ± 0.05 ±0.11 is not large [I], the 
contributions from the flavor and spin violation corrections maybe smear the discrepancy. 

We also explore the possible identification of the D(2760) and D(2750) as the same 3 _ 
state with n = 1, i.e. they are the D% state, the ratio Y(D^D*+n-) = fr° m the heavy 
quark effective theory in the leading order approximation is too large compared with the 

experimental data fr§l§f^h z^R^pi = 0-42 ± 0.05 ±0.11 [1], which again requires the 
flavor and spin violation corrections depressed by the inverse heavy quark mass 1 / uiq are 
as large as the leading order contributions and have opposite sign, such an identification is 
disfavored. On the other hand, the helicity distribution disfavors identifying the -D(2750) 
as the 3" state [lj. We can tentatively identify the (£>(2750), L>(2760)) as the doublet 
(2~, 3~)s with n = 1. 

2 

In this article, we also present the widths for the K and D^rj decays, where 
the strong coupling constants are retained, the predictions can be confronted with the 
experiential data in the future at the BESIII, KEK-B, RHIC, PANDA and LHCb. 



2 We take the approximation r D ( 2 76o) = 1^0(2750)- 
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nLse J F 


Mass [MeV] 


Decay channels 


Width [GeV] 


D* 


1P|2+ 


2460.1 


D* + 7T- 


D+TT- 


0.0543879$; 0.124928$ 


D' 


2s|o- 


? 2539.4 


D*+ir- 


D *o n o 


1.13557$; 0.583137$ 


D*' 


2S$V 


7 2608.7 


D*+tt- 
D*+K- 

D*°TT° 

D*° V 


D+it- 
D+K- 

D°7T 


0.66068$; 0.54317$ 
0.000518592$; 0.106459$ 

0.336747$; 0.276487$ 
0.00841286$; 0.029364$ 


D\ 




7 2763.3 


D*+TT- 

D*+K- 

D*°TT° 

D*°rj 


D+ir~ 
D+K- 

D°7T° 

D° V 


0.339606$; 5.19392$ 
0.0632191$; 1.86912$ 

0.173223$; 2.65247$ 
0.0226441$; 0.508904$ 


D 2 


1£»§2- 


7 2752.4 


D*+TT- 

D*+K- 

D*°7T° 

D*°f] 


D + TT- 

D+K- 

D°7T 

D°r] 


1.27691$; 
0.180643$; 
0.653307$; 
0.069308$; 


D> 2 * 


1D§2" 


7 2752.4 


D*+7T- 

D*+K~ 

D*°TT° 

D*°r] 


D+K- 

D°7T° 
D°7] 


0.221226$; 
0.00413833$; 

0.114719$; 
0.0027123$; 


D 3 


1£>|3- 


7 2763.3 


D*+TT- 

D*+K~ 

D*°7T° 


D+it- 
D+K- 

D°7T° 

D°r] 


0.0907266$; 0.176388$ 
0.00218128$; 0.018115$ 
0.0468994$; 0.0912646$ 
0.00124089$; 0.00618076$ 



Table 3: The strong decay widths of the newly observed charmed mesons with possible 
identifications. 
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nLs(J 



Mass [MeV] 



Ratio 



D* 



1P|2 J 



2460.1 



r(D+n-y 

r(D*+n~) 

r(D+K-) _ 
r(D*+7r-) 



-n M . r(DW) _ n1g . 

= 0.044; r ^; ;2 } = 0.013; 
= 0.001 

= 7.81; r(D*+7r-) = ^-50; 

— f) 51 • r ( D s + A "~) _ n 1 q. 

= 0.067 



25^1- 



7 2608.7 



= 0.82 
0.16; 



r(D* u 7r u ) 
r(D*+7r-) 
r(p°T?) 
r(D*+7r-) " 
r(p;+x-) 
r(D*+7r-) 



7 2763.3 



r(D+7r-) 
r(D*+7r-) 
r(D q) 
r(D*+vr-) 



r(D*"7r") 

r(D*+7r-) 



lo-^y, r(D*+7r-) 
i 50 . r(D*'V») 

r(g*°r?) 
r(D*+?r-) 



0.14; r ^g;'g , = 0.054 



r>9 



1D§2- 



7 2752.4 



0.51; 



r(D*+7r-) 



r(D*+7i-) 



n rr). r(D s * + gfl 

u - oz ' r(D*+7r-) 



D'* 



12?§2" 



7 2752.4 



r(D*+n~) 
T(D+K-) _ 
r(D*+7r-) 



r(g>V>) 

> r(D*+7r-) 

r(p°g) 
r(D*+7r-) r 



0.019; r ^;^ } = 0.012 



= 1-01; f&p) = 0-52; 

= 0.068; rffi+o = °-° 24; 
= 0.014 



£>3 



1D§3- 



Td|T 

1D§2" 



7 2763.3 



= 1.94 
0.20; 



r(D* 



£>2 



1Z)§3- 



7 2763.3 
7 2752.4 



r(D*-j>D+7r- 

r(D 2 -^-P*+7r- 



) _ 



4.07 



^2 

Da 



7 2752.4 
7 2763.3 



r(gg^£>+7r-) _ n 8n 

T(D' 2 *^D*+n-) ~ U ' OU 



Table 4: The ratios of the strong decay widths of the newly observed charmed mesons 
with possible identifications. 



Babar 


CLEO 


CLEO 


ARGUS 


ZEUS 


This work 


1.47 ±0.03 ±0.16 


2.2 ±0.7 ±0.6 


2.3 ±0.8 


3.0 ± 1.1 ± 1.5 


2.8 ±0.8+^ 


2.30 



'-►D+7 



Table 5: The ratio of r ^ D ?( 2 460)°->D*+ -) ^ rom experimental data compared with the 
prediction from the leading order heavy quark effective theory. 
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4 Conclusion 

In this article, we study the strong decays of the newly observed charmed mesons D(2550), 
-D(2600), D(2750) and D(2760) with the heavy quark effective theory in the leading order 
approximation, and tentatively identify the (13(2550), Z)(2600)) as the doublet (0 _ ,1~) 
with n = 2 and (Z?(2750), Z?(2760)) as the doublet (2~,3 - ) with n = 1, respectively. The 
identification of the D(2750) and D(2760) as the same particle with J p = 3~ is disfavored. 
The other predictions can be confronted with the experimental data in the future at the 
BESIII, KEK-B, RHIC, PANDA and LHCb. 

Acknowledgment 

This work is supported by National Natural Science Foundation of China, Grant Numbers 
10775051, 11075053, and Program for New Century Excellent Talents in University, Grant 
Number NCET-07-0282, and the Fundamental Research Funds for the Central Universi- 
ties. 

References 

[1] P. del Amo Sanchez et al, llarXiv: 1009.20761 

[2] A. V. Manohar and M. B. Wise, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 10 
(2000) 1. 

[3] M. Neubert, Phys. Rept. 245 (1994) 259. 

[4] K. Nakamura et al, J. Phys. G37 (2010) 075021. 

[5] S. Godfrey and N. Isgur, Phys. Rev. D32 (1985) 189. 

[6] T. Matsuki, T. Morii and K. Sudoh, Prog. Theor. Phys. 117 (2007) 1077. 

[7] M. Di Pierro and E. Eichten, Phys. Rev. D64 (2001) 114004. 

[8] D. Ebert, R. N. Faustov and V. O. Galkin, Eur. Phys. J. C66 (2010) 197. 

[9] Z. F. Sun, J. S. Yu, X. Liu and T. Matsuki. larXiv:1008.3120l 

[10] X. H. Zhong. farXTv:1009.03591 

[11] P. Colangelo and F. De Fazio, Phys. Rev. D81 (2010) 094001. 

[12] P. Colangelo, F. De Fazio, S. Nicotri and M. Rizzi, Phys. Rev. D77 (2008) 014012. 

[13] P. Colangelo, F. De Fazio and S. Nicotri, Phys. Lett. B642 (2006) 48. 

[14] P. Colangelo, F. De Fazio and R. Ferrandes, Phys. Lett. B634 (2006) 235. 

[15] A. F. Falk, Nucl. Phys. B378 (1992) 79. 

[16] M. B. Wise, Phys. Rev. D45 (1992) 2188. 



9 



[17] G. Burdman and J. F. Donoghue, Phys. Lett. B280 (1992) 287. 

[18] P. Cho, Phys. Lett. B285 (1992) 145. 

[19] T. M. Yan, H. Y. Cheng, C. Y. Cheung, G. L. Lin, Y. C. Lin and H. L. Yu, Phys. 
Rev. D46 (1992) 1148. 

[20] R. Casalbuoni, A. Deandrea, N. Di Bartolomeo, F. Feruglio, R. Gatto and G. Nardulli, 
Phys. Rept. 281 (1997) 145. 

[21] P. Avery et al, Phys. Lett. B331 (1994) 236. 

[22] P. Avery et al, Phys. Rev. D41 (1990) 774. 

[23] H. Albrecht et al, Phys. Lett. B232 (1989) 398. 

[24] S. Chekanov et al, Eur. Phys. J. C60 (2009) 25. 

[25] A. F. Falk and T. Mehen, Phys. Rev. D53 (1996) 231. 

[26] C. G. Boyd and B. Grinstein, Nucl. Phys. B442 (1995) 205. 

[27] Z. G. Wang, Nucl. Phys. A796 (2007) 61. 

[28] S. Ahmed et al, Phys. Rev. Lett. 87 (2001) 251801. 

[29] A. Anastassov et al, Phys. Rev. D65 (2002) 032003. 

[30] Z. G. Wang and S. L. Wan, Phys. Rev. D74 (2006) 014017. 

[31] P. Z. Huang, L. Zhang and S. L. Zhu, Phys. Rev. D81 (2010) 094025. 



10 



